The Algebraic Approach to the Phase Problem for 

Neutron Scattering 

A. Cervellino" and S. Ciccariello^ 

O 

I Consiglio Nazionale delle Ricerche, Istituto di CristallograGa (IC-CNR) 

Campus Universitario c/o Dip. to Geomineralogico, Via Orabona 4, 1-70125 Bari, Italy. 
E-mail antonio . cervellino@ic . cnr . it; Phone +39 080 5442624; 
Fax +39 080 5442591 

, Univcrsita di Padova, Dipartimcnto di Fisica "G. Galilei" & Scz. INFM 

• Via Marzolo 8, 1-35131 Padova, Italy. 

E-mail salvino.ciccariello@pd.infn.it; Phone +39 049 8277173 ; 
Fax +39 049 8277102 



> 

O 



February 7, 2008 



O , 

■rj" , Abstract 

o 

' The algebraic approach to the phase problem for the case of X-ray scat- 

pi |. tering from an ideal crystal is extended to the case of the neutron scattering, 

overcoming the difficulty related to the non-positivity of the scattering den- 
sity. In this way, it is proven that the atomicity is the crucial assumption 
while the positiveness of the scattering density only affects the method for 
searching the basic sets of reflections. We also report the algebraic expres- 
sion of the determinants of the Karle-Hauptman matrices generated by the 
^ ' basic sets with the most elongated shape along one of the reciprocal crystal- 

^ . lographic axes. 
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1 Introduction 



The main crystallographic problem, namely: to determine the electron density of an 
ideal crystal with known chemical composition from its X-ray diffraction pattern, 
is brought to its essence when the atoms are assumed to be point-like because the 
positions of the atoms present in the unit cell are the only unknown quantities to 
be determined. The unknowns' number being finite, it appears reasonable that 
the knowledge of the peak intensities relevant to a sufficiently large portion of 
the reciprocal lattice is sufficient to determine the atomic positions. In fact, Ott 
(1927) and Avrami (1938) first showed that the atomic positions are the roots 
of a set of polynomial equations determined by an appropriate set of reflection 
intensities. This method of inversion of scattering data is known as the algebraic 
approach to the phase problem [Buerger (1960), Hauptman (1991)]. Actually, the 
correct formulation of the algebraic approach is slightly more involved [Navaza & 
Silva (1979), Silva & Navaza (1981), Rothbauer (1994)] for two reasons. Firstly, the 
unimodular roots of the system of polynomial equations, referred to in the following 
as resolvent system, are the positions of the peaks of the infinitely resolved Patterson 
map [Patterson (1939)]. Secondly, for the general case where some of the aforesaid 
peaks have the same projections along one of the three crystallographic axes, the 
resolvent system has to be determined by a more involved procedure than Avrami's. 
These points have been fully clarified in two recent papers [Cervellino & Ciccariello, 
(1996) and (2001)], referred to as I and II in the following. These papers showed the 
existence of many resolvent systems. In particular, the determination of resolvent 
systems is made possible by the positivity of the scattering density, ensured by 
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the fact that we are deahng with X-ray scattering. Very briefly, according to the 
basic paper by Goedkoop (1950), the positivity condition allows us to associate 
to each point of the reciprocal space lattice Zi^ a vector of a finite-dimensional 
Hilbert space 1-L{M). The scalar products of these vectors reproduce the intensities 
of the full diffraction pattern. Each resolvent system is determined by a basic set 
of reflections denoted by B{J\f), i.e. by a simply connected set of J\f reflections such 
that the associated vectors are linearly independent and form a basis of 7Y(A/'). 
The coefficients of the polynomial equations of the resolvent system require the 
knowledge of the peak intensities relevant to all the reflections obtained as difference 
of any two reflections of B{N'). Hence, it is important to select B{Sf) in such a 
way that it is centred on (the origin of reciprocal space) and that its points lie 
as close as possible to the origin. Only when the limiting sphere is large enough 
to contain one of these sets, a resolvent system is known and, after solving it, the 
atomic positions can be determined. The procedure to be followed in order to select 
a basic minimal set of reflections was reported in ref. II, where it was also shown 
how to convert each resolvent system of polynomial equations in three variables 
into a resolvent system of polynomial equations in a single variable. 

The algebraic approach has been successfully applied to solve the structure of 
some real crystals [Fischer & Pilz (1997) and Pilz & Fischer (2000)] and it can be 
implemented to account for experimental errors on reflex intensities (Cervellino & 
Ciccariello, 1999). As a matter of fact, its practical usefulness is severely limited 
by the fact that the degree of the polynomial equations sharply increases with the 
number of the atoms present in the unit cell (Hauptman, 1991). On a theoretical 
ground, the approach looks however quite interesting for its rigorous conclusions 
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and its far reaching implications since the algebraic approach is intimately related 
to other classical issues of mathematical-physics (see the introductory section of 
II). 

The aim of this paper is to report on the extension of the algebraic approach 
to the case of neutron scattering. As already mentioned, the presently known 
formulations exploit the positiveness of the scattering density, a condition generally 
not fulfilled in the case of neutron scattering due to the fact that some atomic 
species have negative scattering lengths. On this basis, one rightly wonders whether 
the mentioned results - in particular the property that the full diffraction pattern 
can be reconstructed from the knowledge of the intensities relevant to a finite 
set of reflections, i.e. the "difference" set generated by a basic set of reflections 
- do apply to neutron scattering or not. We shall show that the answer to this 
question is affirmative^. In order to prove this statement, it is necessary to relax 
the positiveness condition. Hence, the plan of the paper is as follows. In §2 we 
report the basic equations of the algebraic approach and the finite vectorial space 
Ti.{Af) will be introduced on the basis of simple quantum mechanical notions. Based 
on the results proven in Appendices A and B, in §3 we generalize the algebraic 

approach to neutron scattering in the case of two-dimensional crystals and in §4 we 

^To the authors' knowledge, Navaza & Navaza (1992) already gave a positive answer to this 
question. However, these authors explicitly recognized to see no rigorous way for demonstrating 
the reconstruction procedure in three dimensions, that is the most interesting point. This difficulty 
is related to that of singling out a basic set of reflections from the set of the observed ones. In 
this paper we overcome this difhculty by showing how the isolation of a basic set of reflections can 
be carried out also in presence of a non-positive scattering density. By so doing, we generalize 
the results obtained in papers I and II and based on the positeveness of the scattering density. 
It should also be remarked that the aforesaid proof of the reconstruction property requires no 
probabilistic assumption. Therefore, our conclusions are more general than those obatined by 
Hauptman (1976) with the probabilistic approach that, for practical appplications, is by far the 
most useful one (see, e.g., Hauptman and Langs (2003)). 
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sketch the generahzation to the three-dimensional case and report our conclusions. 
Appendix A illustrates a new procedure, not requiring the positivity assumption, 
for singling out a principal basic set of reflections, where principal means that the 
basic set has the most elongated shape along one axis of reciprocal space. Appendix 
B deals with the derivation of the algebraic expression of the determinant of the 
Karle-Hauptman matrix associated to a principal basic set of reflection. 

2 Basic results of the algebraic approach 

The formulation of the algebraic approach, reported in I and II, assumed positivity. 
We will now retrace our steps through the theory in order to make the necessary 
changes to allow for non-positive scattering densities, as it happens with neutrons. 
We continue to assume that the unit cell contains point-like atoms. Its scattering 
density has the following expression 

N 

Pceiiir) = J2Zj6{r-rj) (1) 

where is the three-dimensional (3D) Dirac function, tj the position of the jth 
atom and Zj the atomic number or the scattering-length of the jth atom, depending 
on whether one considers X-ray or neutron scattering. The two cases differ because 
the 2j^s are positive integers in the case of X-rays and only real numbers in the 
case of neutrons^. Thus, the positiveness of the scattering density is generally not 

ensured in the neutron case. But Iobs,h - the intensity observed at reflection h - is 

^It is understood that absorption and other experimental effects are either absent or corrected 
for. 
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in both cases the square modulus of the Fourier transform of (1), i.e. 

2 N 



N 



E 44e^^-'^-(^^-^). (2) 



Each vector (r^ — r^) can be brought within the unit cell by adding to it a vector 
nij with components equal to or -1, so as to write 

Tj - Ffc + nij- fc = S. (3) 

As (j, A;) runs over its A^(A^ — 1) values, we label the different 5's, defined by (3), 
by j and we denote by JV' the number of the different 5j's. Moreover, we denote by 
£j the set of pairs (j, k) such that (r^ — Vk) defines the same Sj after applying (3). 
Then, the second sum on the right hand side (rhs) of (2) becomes 

^g^2.h.?, ^ 2jZ, (4) 



After setting 



and 



E (5) 



(i,fc)G£j 



Eq.(2) reads 



N 

i=i 

/h = E^^ie^^^'^-'X (7) 

where A/" is the number of the i^j's different from zero. [In the case of neutron 
scattering. A/" can be smaller than Af' because the negativeness of some Zj^s can 
make some i^j's equal to zero.] Eq.(7) shows that the Jh's, the "subtracted" peak 
intensities defined by Eq. (6), are the Fourier transforms of the scattering density 
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relevant to the Patterson map 

P.M = T.^sK^-^5)- (8) 

i=i 

This consists of M scattering centres located at 5i, . . . , with weights or "charges" 
equal to z/i, . . . , uj^, and the positiveness of the weights is ensured only in the case 
of X-ray scattering. Moreover, Eqs. (6), (7) and (8) make it evident that the 
knowledge of all the observed intensities Iobs,h. only determines quantities 5j's and 
i^'j's, i.e. the scattering density of the Patterson map. Assuming the latter quantities 
known, Eqs (3) and (5) can be inverted to determine the atoms' positions ri, . . . , tat 
by the procedure reported in §3.2 of I. This deconvolution of the Patterson map 
involves a finite number of operations. In this way, all the atomic configurations 
that reproduce the observed diffraction pattern are determined. Hence, the difficult 
problem to be solved is to find out the set of Eq.s (7) that uniquely determine 
5 J and z/j for j = 1, . . . , A/". The solution of this problem requires, firstly, the choice 
of an appropriate set of h values that determine the equations to be solved and, 
secondly, a procedure able to solve the resulting set of non-linear equations. 

For X-ray scattering, the solution of the first problem is achieved by introducing 
the Goedkoop (1950) lattice of vectors, which is a subset of a finite-dimensional 
Hilbert space. Unfortunately, this step requires that all the i^j's are positive and, 
therefore, it cannot be extended to the case of neutrons. However, by using some 
notions of elementary Quantum Mechanics, we show now that in both cases it is 
possible to introduce a finite-dimensional Hilbert space and, within the latter, a 
lattice of vectors in such a way that the scattering density (8) and the "subtracted" 
intensities (7) are two different representations of a single hermitian operator. 
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To this aim we recall that the position and momentum operator, respectively 
denoted by TZ and P, have eigenvectors |r) and |p) whose eigenvalues r and p span 
the full 3D space B?. Consider now the eigenvalues p equal to — 27rh, h being a 
triple of integers, and put |h) = |— 27rh). As h ranges over the 3D lattice Z^, the 
set of |h)'s defines a lattice of vectors lying within the infinite-dimensional Hilbert 
space Ti spanned by the eigenvectors |p) or |r). Introduce now the linear operator 

Q^j:\6s)uM^ (9) 
i=i 

where is the eigenvetor of TZ with eigenvalue (5j equal to the position vector 
of the jth scattering centre. Due to the property (r|r') = 5{v — r'), the matrix 
elements of Q with respect to the eigenvectors of IZ are 

(v\Q\v') = 6{v-v')Y.^~,6{v-6s) (10) 
i=i 

At the same time, the matrix elements of Q with respect to the lattice vectors |h) 
are 

A?" 

(h|Q|h') = (27r)-3^z/je'2-^^('^-'^'\ (11) 
j=i 

where we used the property that (p|r) = e*P '"/(27r)'^/^ and units such that h = 1 
(Messiah, 1959). Comparison of (10) with (8) shows that the scattering density 
(8) coincides with the diagonal matrix elements of Q (leaving aside the divergent 
factor related to the value 5(0) of the 1st Dirac function). On the other hand, the 
comparison of (11) with (7) shows that all the subtracted intensities (7) are {2ttY 
times the matrix elements of Q with respect to the lattice vectors |h). Moreover, 
Eq.(9) shows that the "charge density" operator Q is determined only by the M 
eigenvectors . . . , \5fj) of IZ with eigenvalues equal to the positions of the N 
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scattering centres, and by the M real numbers z/i, . . . , z/_^ equal to the weights of 
the scattering centres. Hence, we can restrict ourselves to the finite-dimensional 
Hilbert space T-C{J\f) spanned by the vectors \6i), . . . , \6j^) and defined as 

r ^ - 1 
^ i=i ^ 

Vectors \6i) , . . . , \6j^) obey the orthonormality condition^ 

0j'\^j)=Sj',h i,i' = i,...,Ar, (12) 

being the Kronecker symbol, as well as the completeness relation 

J:\63){63\ = 1. (13) 
i=i 

In order to preserve the validity of (11), we still need to assume that 7i(A/') contains 
a lattice of vectors |h) [not to be confused with |h) or with the eigenvectors of V, 
see the following Eq. (17)]. To his aim, it is sufficient to put 

A/" 

|h) = ^ e-'2-h-5i|^^)^ Vh G (14) 
i=i 

After taking the scalar product with {6j' \ one gets 

|h) = e-*2-h.5,, ^ Vh G Z^, j' = 1, 2, . . . , AT. (15) 

/^From the above two relations it follows that vectors |h) are no longer orthogonal 
since from (14) and (12) one gets 

A?" 

(h|h') = ^e^^--^^-^^"^') = (h + m|h' + m), ^h,h',meZ\ (16) 
i=i 



^By so doing, the previous normalization {Sf\6j) = 6{5j' — 5j) has been scaled to {Sf\6j) — Sj'j. 
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with (h|h) = M . This property is not surprising if one observes that |h) and |h) 
are related as follows 

|h) = (27r)^/^f:|5;)(d;-|h) (17) 
i=i 

so that |h) is the projection of |h)(G Ti) into 1-L{M) and, therefore, it is no longer 
an eigenvector of V. Now it is important to note that 

i=i 

Thus, on the one hand, all the matrix elements of Q with respect to the lattice of 
vectors |h) reproduce the full diffraction pattern. On the other hand, the diagonal 
matrix elements of Q with respect to the basis vectors are the weights of the 
scattering density (8). In this way, it has been shown that: i) both for X-ray and 
for neutron scattering it can be introduced a finite-dimensional Hilbert space 7i(A/') 
spanned by the M eigenvectors of TZ associated to the position vectors of the M 
scattering centres, ii) within 1-L{M) it exists a lattice of vectors = {\\\) \ h. E Z^} 
with |h) defined by Eq. (14), iii) it exists a hermitian linear operator Q whose matrix 
elements with respect to the basis vectors and to the vectors of the vectorial 
lattice Z^ yield all the weights of the scattering density and all the subtracted 
intensities /h, respectively. 

We are now left with the problem of determining A/", the |5j)'s and the i^j's 
knowning an appropriate number of Jh values. Before tackling with this problem 
in the following sections, we report some interesting consequences of the aforesaid 
vectorial structure underlying the phase problem for an ideal crystal. The first, 
related to Eq. (16), shows that (h'|h) = (h' + m|h + m), Vh, h', m G Z^ . Thus, 
the scalar product of any two vectors of Z^ does not change if the associated lattice 
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points are translated by an arbitray vector m E Z^. The same property applies to 
(h'|Q|h) because by Eq. (18) one gets 

/h'-h= (h'|Q|h) = (h' + m|Q|h + m) = (h'-h|Q|0), V h, h', m G (19) 

Second, cannot contain more than Af linearly independent vectors because it 
is a subset of 7i(A/'). In Appendix A it will be shown that exactly contains 
A/" linearly independent vectors. Thus, if we denote one set of these vectors by 
|ki), |k2), . . . , |k_^), we can write 

AT 

|h) = E^h,,m-), VheZl (20) 
i=i 

Taking the scalar product of the adjoint of this equation with vector |0) and using 
Eq. (19) one obtains 

AT _ Af _ 

/h = E^hj(k,IQ|0) = E^h,,^k, Vh e Z^ (21) 
j=i i=i 

where the overbar denotes the complex conjugate. Eq. (21) shows that any sub- 
tracted intensity Jh is a linear combination of the J\f intensities Jk^, . . . , /k^. The 
matrix elemets obey to a set of relations. The first is obtained by taking the 
scalar product of (20) with {6j\ and reads 

^-i2nh-5, ^ ^ ^^^^g-i2.k,.5,^ yj^ ^ ^3^ j = 1 , 2, ■ ■ ■ , AT. (22) 

e=i 

After introducing an A/" x A/" matrix (V) with Vj/ = e~'^'^'^^''^3 ^ the previous equations 
becomes 

EVj,,Ah,, = e-^2-'^-^^i, (23) 
which, as it will be shown later, can formally be solved as 

Ah/ = E(V-^),je-^2-'^-^^ ^24) 
i=i 
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The second is obtained by substituting Eq. (20) in Eq. (19), using Eq. (21) and the 

hnear independence of |ki), . . . , One finds 

M M _ _ 

^-h'i = E E ^h,i^h'j'^,-v,z Vh, h' 1 = 1,2,. (25) 

3 Neutron scattering 

Given their lengthy nature, we report in Appendix A the proof of the existence of 
the principal basic set of vectors |h) along a crystallographic axis, and in Appendix 
B the analytic expression of the determinant of the Karle-Hauptman matrix asso- 
ciated to this principal basic set of vectors. It is stressed that these results, worked 
out in the two-dimensional case {D = 2), apply both for X-ray and for neutron scat- 
tering, because we never require the posivity in proving them. Moreover, they easily 
generalize to the three-dimensional case. In this section we show how to single out 
a basic set of refiections in the case of neutron scattering, knowing an appropriate 
number of subtracted intensities Jh. [For simplicity we shall still confine ourselves 
to the case D = 2.] Firstly, we recall that the results, obtained in the previous 
section as well as in Appendices A and B independently of the positeveness, show 
that: i) it is possible to introduce a Hilbert space 7^ (A/") having as an orthonormal 
basis the set of the (scaled) eigenvectors \6j) of 71 with eigenvalues equal to the 
position vectors 6j of the A/" scattering centres of the (infinitely resolved) Patterson 
map determined by the full diffraction pattern, ii) within 7i(A/') it is possible to 
extract a lattice of vectors |h) defined by (14), iii) the (subtracted) intensities 
defined by (6) are the matrix elements of the charge density operator Q [defined by 
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(9)] with respect to vectors of Z^, while the charges z/j are the eigenvalues of the 
eigenvectors of Q, iv) the knowledge of the eigenvalues (5j allows us to determine 
the subset X [defined by (39)] of Z"^ and the translation of 1 by (—1, —1) yields the 
principal basic set i3(a*) of Ti^M) along the reciprocal crystallographic direction 
a* and, finally, v) the algebraic expression of the determinant of the KH matrix 
associated to i3(a*) is given by Eq. (69). Even though results i)-iv) were essentially 
obtained in papers I and II, it is stressed that they are now extended to the case 
of neutrons. Moreover, the above presentation makes the introduction of T-C{Af) 
clearer and shows that the i^j's and the Jh's are the matrix elements of a single 
operator (Q) with respect to two different sets of vectors. As yet, however, the 
aforesaid generalization is practically useless, because we do not know the (5j's and 
JV. Since the only known quantities are the subtracted intensities Jh, the determi- 
nation of the 5j's, z/j's and A/" must be carried out in terms of the Jh's. Hence, the 
search of a basic set must be performed in terms of these quantities. In papers I and 
II, we reported the procedures for carrying through such a search. Unfortunately, 
they only apply to the case of X-ray scattering because they exploit the positiveness 
of the charge density operator Q, a condition fulfilled only in the case of X-rays. 
Very briefly, as shown in I and II, the simplest search of a basic set proceeds as 
follows. One starts from the set of reflections B2 = {(0)0), (1,0)} and one eval- 
uates the determinant of the associated KH matrix (^P[i32]^. If det^P[;B2]^ 7^ 0, 
one "enlarges" B2 by "adding" to it the next reflection (2,0) so as to have the new 
set B3 = {(0, 0), (1, 0), (2, 0)}. Then one evaluates det^P[i33]^ and if this value 
is different from zero one adds the next reflection (3,0) to so to have the en- 
larged set B4. As far as the vectors associated to the reflections of the considered 
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set are linearly independent, the positiveness of Q ensures that the corresponding 

KH matrix has a strictly positive determinant. On the contrary, if the vectors are 

linearly dependent the determinant is equal to zero. This property is easily shown 

as follows. Assume that Bm. = {hi,h2, . . . ,hm} and consider the m x m matrix 

(Qm) having its (r, s) element given by 

AT ^ ^ AT _ 

(hr|Q|h^) = /h,-h, = Yl {^r\Sj)i^j{Sj\hs) = Vj,ri^jVj,s, r, s = 1, . . . , m (26) 

i=i j=i 

where Eq. (15) and the definition of Vj,/ reported below Eq. (22) have been used. 
Eqs.(26) can be written in matricial form as 

[Qm) = (Vt) . (.) ■ (v) (27) 

where ^V^ is now an J\f x m matrix (with m < J\f), (y'^^ its hermitan conjugate 
and (v) an M X M diagonal matrix with its {]■,]') element equal to i^jSjj. The 
determinant of matrix (Qm), evaluated by Bezout's theorem (Gantmacher, 1966), 
yields 

det(Q^) = Y 

i<h<h<--<Jm<Af 

where (Vjj....j^) denotes the mxm minor formed with the jith,...,jmth row of (V). 
In the case of X-rays the non-negativeness of the addends present in the above sum 
implies that det(Qm) 7^ unless all the quantities det(Vji,...j„) are equal to zero 
for 1 < ji < 72 <■■■< im < A/". The latter conditions are verified if and only 
if the considered vectors |hi), . . . , |hm) are linearly dependent. In fact, when this 
condition is fulfilled, the generic vector 

m 
r=l 
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det(Vj„ 



(2J 



can be equal to zero with some a^-'s different from zero. Then, the scalar products 
with {6j\ yields 

m m 
= H Oir{Sj\hr) = arVj,r = 0, j=l,...,J\f. 

r=l r=l 

This being a set of Af linear homogeneous equations in the unknown ai, . . . , am, 
the Rouche-Capelli theorem ensures that a non-trivial solution exists if and only if 
the rank of the A/" x m matrix (V) is smaller than m, i.e. if the determinants of all 
the m X m matrices contained in (V) are equal to zero. In this case, one finds that 
det(Qm) = and the vectors of the considered set Bm are linearly dependent. On 
the contrary, if one of the mxm matrix contained in (V) is non singular, all the a^'s 
are equal to zero, the vectors |hi), . . . , \hm) are linearly independent and det(Qm) 7^ 
0. Coming back to the search procedure of a basic set, the aforesaid property makes 
it clear that the enlargement procedure comes to a halt when det(V[Bfn]) = 0, i.e. 
when one finds a "KH zero". This condition must certainly occurr because the 
number of the linearly independent vectors cannot exceed A/". Hence, according to 
the analysis reported in Appendix A, we find that rh = fii = M and the vectors 
associated to the set B^-^ = {(0, 0), ... , (M — 1, 0)} are linearly independent. Then, 
we enlarge the previous set by adding to it, step by step, the reflections lying on 
the next upper row starting with (0, 1). The next KH zero is found when we "add" 
the reflection (/i2, 1). We move now to the next upper row and we start by adding 
the reflection (0,2) to the set B^,,,^^ = {(0, 0), . . . , (/^i - 1, 0), (0, 1), . . . , (/i2 - 1, !)}• 
The search of the basis set is accomplished when we arrive at the point where the 
inclusion of the reflection (0, mi) leads to a KH zero. By construction, the resulting 
basis set is simply connected. 
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In the case of neutron scattering, this analysis is no longer possible because 
the finding of a KH zero {i.e. det(Qm) = 0), during the enlargement procedure, 
does not allow us to infer that the associated vectors are linearly dependent. This 
appears evident from Eq. (28): here, each factor related to the charge product is 
not ensured to be positive so that the condition det(Qm) = does not imply that 

2 

all the quantities det(Vji,...j-„) equal to zero. Thus, in order to extend the 
search procedure of a basic set to the case of neutron scattering, we must introduce 
a positive definite operator whose matrix elements with respect to the vectorial 
lattice are known in terms of the observed scaled intensities I^. To this aim, 
denote by Sobs the set of the observed reflections and denote by Si the largest subset 
of Sobs such that for any two reflections and h^' in Si it results (hj, — h^') G Sobs- 
We denote by A/i the number of reflections contained in Si and we assume flrst 
that Si is large enough to contain at least one basic set so that A/i > A/". Consider 
now the linear operator^ 

Qs, = QE|hr)(h.|Q. (29) 

r=l 

This operator is hermitian and positive deflnite. The flrst property is evident. 
To show the second, we consider the expectation value of with respect to an 
arbitrary vector \a) G Ti^N). One flnds that 

Ml , 



Qs,\a) = ^ (a|Q|h,) 



r-l 



This expectation value can be equal to zero either if Q\a) is perpendicular to all 

the |hr)'s for r = 1, . . . , A/i or if Q\a) = 0. The flrst condition is impossible unless 

''The introduction of this quantity is suggested by the procedure followed by Silva & Navaza 
(1981) and Navaza & Navaza (1992). 
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Q\a) = because Si is assumed to contain a basic set of vectors. We are left with 
the condition Q\a) = 0. This imphes that \a) is eigenvector of Q with eigenvalue 0. 
But this condition is impossible because the eigenvalues of Q are all different from 
zero as it appears evident from Eq. (9). Hence, {a\Qs^\a) > V|a) ^ and the 
positivity of is proven. The matrix elements of with respect to the vectors 
of the lattice are 

(h|Q5jk) = ^/h-h,./h.-k. 

r=l 

If h, k G iSi, the matrix elements of are fully known and will be denoted as 

Ml 

Jh,,h„, = (hi|Q5i|hr„) = Zl^h,-h./h.-h„, l,m = l,...,Afi. (30) 

r=l 

[It is noted that the i7h,,h„'s are symmetric since they obey the relation i7h,,h„ = 
i/h,„,h; that follows from the Friedel property valid for the subtracted intensities, 
i.e. Jh = I-h-] At this point, the search of a basic set becomes possible acting as 
follows. We start from the set of reflections B2 = {(0,0), (0, 1)} and we evaluate 
the determinant of the matrix whose elements are the matrix elements of Qs^ 
between the vectors associated to 82- These matrix elements are known owing 
to (30). For simplicity, this matrix also will be called a Karle-Hauptman matrix, 
even though its matrix elements are the Xh^^h^'s instead of the /h;,h,„'s- If the 
determinant of this KH matrix is different from zero, we enlarge B2 by adding 
to it a further reflection chosen either by the procedure described above or by 
one of the other procedures reported in II. Assume that Bm is the first set, found 
during the enlargement procedure, such that the determinant of the KH matrix 
(with elements Ihi,h^) is equal to zero. Since Qs-^ is a positive definite operator, 
the vectors |hi),..,|hm) associated to the reflections of Bm are linearly dependent. 
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Moreover, the vectors |h) considered in this paper do not depend on the charges z/j 
[see Eq. (14)]. Actually, as explained in footnote 6, they refer to positive charges. 
Thus, property 1 proven in sect. 4 of II applies. Therefore, in the subsequent 
enlargement procedure, we must discard all the reflections relevant to the quadrant 
defined in property 1. The search of a basic set will be accomplished once the 
resulting set cannot further be enlarged by the adopted procedure of enlargement. 
If this happens, one concludes that the considered Si set and, consequently, the 
underlying limiting sphere are large enough to contain a basic set. Before discussing 
the consequences of this result, we need to say what happens when Si is not large 
enough to contain a basic set of vectors. In this case, we can always denote by 
|hi), . . . , IIim) the M(< A/i) linearly independent vectors contained in Si that are 
closer to the origin of reciprocal space. The remaining vectors |hA/+i), . . . , \h/^^) of 
Si take the form 

M 

Ihr-) = E«ni|hi), r = (M + 1),...,M. 

r=l 

The aforesaid linearly independent vectors can be singled out by the centred square 
procedure reported in §5 of II because Qs^ is positive definite in the subspace 
spanned by |hi), . . . , Ih^/). By so doing we shall find some KH zeros, but the KH 
matrix (with elements Xh;,h„) relevant to the set of vectors Bm = {|hi), . . . , jh^)} 
certainly is non-singular. On the contrary, the KH matrices relevant to the (A/i— M) 
sets of vectors Bm U \hr), with r = (M + 1), . . . , A/i, are singular and all the h^, 
with r = {M + 1), . . . , A/i are KH zeros. However, the locations of these zeros must 
be such that the configuration of Bm is not that of a basic set, in the sense that 
the locations of the KH zeros is such that an "enlargement" of Bm by the centred 
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square procedure is possible whenever one could dispose of a set S[ D Si (as it 
would happen with a larger limiting sphere). One concludes that the configuration 
of the KH zeros associated to Bm is alike to that found in the previous case, i.e. 
when Si contains A/" linearly independent vectors and a basic set is not found. 
Hence, only two cases are possible: either Si is large enough to contain a basic set 
or it is not. In the second case, the observed diffraction pattern does not allow to 
solve the phase problem. In the first case, it does. In fact, as already stressed, the 
singled-out basic set is also a basic set for the subtracted intensities Jh and we can 
use all the results found in papers I and II for the case of X-ray scattering^. In 
particular, the knowledge of the basic set and of the further refiections where we 
have found a KH zero allow us to determine the associated complete set (C) and 
complementary set (JF) of refiections. The first consists of the refiections differences 
of any two reflections of the basic set B/j-, and the second of the reflections (not 
contained in the complete set) that are differences of the reflections associated to 
the KH zeros with the reflections of the basic set. The recursive determination of 
the intensities relevant to the reflections not contained in C U JF proceeds along the 
lines described in I and one flnds that also in the case of neutron scattering the 
full diffraction pattern is determined by the knowledge of the intensities associated 
to the flnite set of reflections C U JF. Moreover, each KH zero states that the 
added vector to Bm is a linear combination of the vectors of Bm- The coefficients 

of the linear combination are obtained by solving the associated linear equations 

^Strictly speaking, the determination of matrix (T) as reported in Appendix D of I is not 
possible in the case of neutron scattering. Hence, the matrix (7^), instead of being evahiated by 
Eq. (1.2.41), must be evaluated by the first equality in Eq. (1.2.45), i.e. by inverting the KH 
matrix (with elements /h,,h„) associated to the basic set of reflections. 
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whose matrix of coefficients is the KH matrix associated to Bm- The scalar product 
with {6j\ yields a polynomial equations in the variables e"*^'^^^ and e~*^'^^^ (we 
recall that we are considering the case D = 2). The system of these polynomial 
equations has 2A/' common unimodular roots that are related to the positions of 
the J\f scattering centres as it is just specified. It is interesting to note that, by 
following the procedure reported in sect. 6 of II, the aforesaid system of polynomial 
equations in two variables can be converted into a set of polynomial equations in a 
single variable, so as to make the solution of the problem simpler. 

4 Conclusion 

The results reported in §3 and appendices A and B refer to the two-dimensional 
case. Their extension to the case D = 3 is rather straightforward by following the 
lines illustrated in paper II. In particular, we write now 6j as (xj, ^j^^^) where 
index t labels the different projections of all the 5j's along a = x, j the different 
projections of the (5j's that have the same x projection along h = y, and £ the 
different projections along c = ^ of the 5j's that have the same x and y projections. 
Then the set X becomes 



l<i<%j, 1<J<P^, l<^<M}, (31) 



with Af = J2fii I]j=i Qi,3- Here label t is assigned in such a way that pi > P2 > 
■ ■ ■ > Pm- Then, for each i, label j is defined in such a way that g^^i > ■ ■ ■ > qi,p^, 
and similarly for i. After putting 

= e-^2""% 7],,, = e-'^'^y-^ and = e"*''^^-^ (32) 
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a generic vector |h) with h = [h, k, I) takes the form 

i=\ j=i £=1 



(33) 



By the same analysis carried out in Appendix A, one finds that the basic set of 
reflections consists of the refiections (/i, fc, /) determined by the conditions: < / < 
— 1), < < (pi+i — 1) and < /i < (M — 1), so that i3(a*) coincides with 
X shifted by (-1,-1,-1). The generic element Vhj of the associated M ^ M matrix 
(V) now reads ii^f]i,j^C,\ ] e. ^^d, generalizing the procedure followed for obtaining 
Eq (66), one finds that its determinant is 

■ M Pt 



det(V) 



n n n ( c.,.. - c.,..) 

1=1 ]=1 l<si<S2<qi,] 



n [in- 

l<'tl<«2<A/ 



M 



n n I "^^2 - Vi,ri 
1=1 l<ri<r2<Pi 



(34) 



where it has been put = Z]j=i Qij for z = 1, . . . , M. Since Eq. (68) holds also true 
in the 3D case, one concludes that the determinant of the KH matrix associated to 
the 3D principal basic set of refiections defined above has the following algebraic 
expression 

2min{P,2,P,J 



det D[S(a*)] 



M Pi 

nn n 

4 = 1 J=l l<Sl<S2<(Jij 

M 



n n 

1=1 l<ri<r2<Pi 



2min((3i,r2igi,r2) 



n 

1<11<22<M 

AT ■ 

n-i 



it2 



i=i 



(35) 



In conclusion, it has been shown that the algebraic approach applies both to X- 
ray and to neutron scattering. In this way it appears clear that the atomicity is 
the crucial assumption, while the positiveness of the scattering density is from a 
theoretical point of view not important. It only makes the search of the basic set 
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faster, because the search can be carried through with the subtracted intensities /h- 
On the contrary, in the case of neutrons, the procedure is shghtly more involved, 
because one must first select the largest Si subset within the observed diffraction 
pattern and by these intensities to evaluate the jTh^.h^'s for all the and of 
Si. Then, the search of a basic set is performed, as in the case of X-ray scattering, 
using the i7hr,hs's. Finally, after finding a basic set, one proceeds with the /h's and 
with the found basic set to reconstruct the full diffraction pattern and to determine 
the atomic positions as in the case of X-ray scattering. 
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A Principal basic sets for the 2-D case 

We show now that contains different sets of M hnearly independent vectors 
and we report the procedure for selecting one of these sets. Actually, this proof 
is immediately achieved by the procedures illustrated in I and II if we assume to 
know the quantities 

/h,k = (h|k) = e^'^^*^-^)-^^^ = Jh-k (36) 

[where the rhs follows from Eq. (16)], because the involved vectors belong to 
and, therefore, belong to a Hilbert space so that the assumptions, made in I and II, 
are obeyed^. We prefer however to proceed differently in order to make it evident 
that the geometry of the scattering centres determines the principal basic set of 
reflections. To this aim, we first observe that all the equations reported in §2 hold 
also true if we restrict ourselves, for greater notational simplicity, to the case of a 
2D space. We remark that, even though the locations of the M scattering centres 
are different from each other, it can happen that the distinct projections of the (5j's 
along one of the crystal axes are only M with M < Sf . Hence, we shall denote the 
distinct projections along x = a by xi, . . . , xm- Consider now those (5j's that have 
x-projections equal to, say, xi. Since these 5j's are different, their projections along 
axis y must differ among themselves and we shall denote their number by nii > 1. 
In this way, the considered 5j's can be written as [(xi, |/i), (xi, 1/2)5 • • • , (a^i, ymj] 
and, in general, we have 

. . .Jsr^ = l^{xr,yr,s) I s = 1,2, . . . ,m^, r = 1, . . . ,M, I (37) 

^The comparison of (36) with (7) shows that the /h-k's can be considered as the scattering 
intensities relevant to a set of J\f scattering centres, located at 5i, . . . , 6j^, with positive charges 
all equal to one. 
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Further, in labeling the different Xr, we choose the label in such a way that rrii > 
777,2 > ■ ■ ■ > f^M- It is observed that points {xr, Urs) can be mapped into a subset 1 
of the Z'^ lattice defined either as 



or as 



J={(r,s) 



1= {(r,s; 



s — 1, . . . , m^, r = 1, . . . , M} (3J 



r = l,...,/is, s = l,...,mi} (39) 



In the first case, is the number of points lying on the rth column, while in the 
second //^ is the number of points lying on the sth row. In both cases, we have 
Erli"^r = TJ^=ilJ^s = A/". Writing |h) as \h,k) with h and k integers, Eq. (14) 
reads 

M rrir 
r=l s=l 

or, after putting 

^e"^^"^-- and r/,,, = e-^^'^"--, (40) 



as 

M nir 

1^'^) = im^r-V,/kr,yr,s)- (41) 
r=l s=l 

Consider now the polynomial 

M M 

PA^) = I[i^-^r)=T.»M,sZ' (42) 
r=l s=0 

with 



-1)''"^ E s = o,...,(m-i), (43) 

l<sx<S2"-<s^j_g<A/ 
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and = 1. Eq.(42) implies that P^,(^r) = for r = 1, . . . , M. From Eq. (42) 

one gets 

M-l M-1 

Cr-"^ = - E «M,.e/ = E Pm,m,s^/, r = 1, . . . , M (44) 

s=0 s=0 

and substituting these relations in |M, k) one obtains 

M-l M rrir M-l 

\M,k) = - E "a/,pEE^^^'^'-:^''I^'-'^^:^) = E pM,M,s\s,k), (45) 

p=0 r=l s=l s=0 

which shows that |M, /c) is a linear combination of vectors \s, k) defined as 

M nir 

\s,k) =J2T.^r"Vr/\Xr,yr,p), S = 0, . . . , (M - 1) . (46) 
r=l p=l 

Moreover, after multiplying both sides of Eq.(44) by and using again Eq.(44), 
one gets 

M-2 M-l 

= - E - «M,(M-.)e.^ = E r = i, . . . , m, (4?) 

s=0 s=0 

with 



a _ I '^Af,(A/-l)'^i\/,0' ^ ^' I'^c 

A^A/,(A/+i),s 1 -Q, Ma a if s = 1 (M-l) ^ 

I "a/,s-i ^ "a/,a/-i"a/,3 ' -^5 • • • ! V^"^^ -^y- 



Eq. (47) shows that also is a linear combination of • • • , ^r*^ ^ with coef- 

ficients Pj^j j^^^^j ^ specified by Eq. (48) and obtained by a recursive application of 
Eq.(44). After dividing the equation P^^(^r) = by ^r~^, one obtains that 



A/-1 A/-1 



^ = -aj,0 E "A/.s+l^r" = E /^A/.-i.s-^r-", r = 1, . . . , M. 

Iterating the procedure one finds that 

M-l 

= E PM,Ur% r = l,...,M, \/k>M, \/k<-l. (49) 

s=0 
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/^From these relations it follows that 

M-l 

\h,k)=Y^ Pmas\s, k), \/h>M,\/h<-l. (50) 
Consider now the polynomials 

rrir 

PmA^) = Iliz-Vr,s), for r = l,...,M. (51) 

s=l 

By the same procedure, one proves that 

rrir — l 

^r/ = Pmr,k,qVr/, T = 1,...,M, s = 1, . . . , m^, V/c > m,., V /c < -1 (52) 

g=0 

where coefficients (3mr,k,q are iteratively determined in terms of the coefficients that 
define the polynomials (51). We put now 

rrir 

\Xr, k) = Tjr^s 'l^ry yr,s) y k = 0, 1, . . . . (53) 
s=l 

^From Eqs (12) and (40) it follows that 

M 

{Xg, k \Xr, k) = 6s^r f]r,s (^4) 
s=l 

and one concludes that vectors \xr, k) and |xs, k') are linearly independent if r 7^ s. 
The same happens for the vectors 0), ■ ■ ■ (m^ — 1)). This property im- 
mediately follows from Eq. (53), because the rhs of the equation involves mr 
linearly independent vectors at fixed r and the matrix of coefficients ?7r,s'^, with 
/c = 0, . . . , (m^. — 1) and s = 1, . . . , m^, is a Vandermonde matrix with determinant 
equal to ni<p<g<mr ('7'-,ij ~ 'Ur,p)^ which is certainly different from zero because the 
T^r^q s are all different among themselves. Combining Eqs (41) and (53) one finds 
that 

M 

\h,k) = Y,^r''\xr,k). (55) 

r=l 

26 



If /c = 0, all the vectors on the rhs are linearly independent as r ranges from 1 to M. 
If we let h range in [0, . . . , (M — 1)], we have M linear relations and the coefficient 



<i<j<M 



matrix is a Vandermonde matrix with determinant equal to 
Thus, the linear independence of 0), . . . , |a;M, 0) ensures the linear indepen- 
dence of the vectors |0, 0), . . . , |(M — 1), 0). If rriM > 1, by the same proce- 



dure and using the linear independence of the sets of vectors 



\xi, k), . . . , \xm, k) 



for k = 0, . . . , {rriM — 1), one finds that each of the following sets of vectors 



|0,£),...,|(M-1),. 



with 



0, . . . , {rriM — 1), is a set of M linearly inde- 



pendent vectors. Moreover, the vectors 



|0,^),...,|(M-1),0 



are linearly inde- 



pendent from the vectors [|0, j), . . . , |(M — 1), j)] if < ^ 7^ j < (m^/ — 1) due to 
Eq. (54). Consider now the mxm Vandermonde matrix I V(m, ^) ) with elements 



(56) 



As far as m < M, the aforesaid matrix is non-singular and endowed of an inverse 
denoted by I V~^(m, ^) ). Then, one can write 



M 



.q) = Y.^-\,s{M,0\is-l),q), l<r<M, < g < (m^, - 1). (57) 



s=l 



Assume now that mM-i > friM and consider those values of k such that ttlm < k < 
rriM-i- Eq- (41) can be written by (52) as 

M—l rrir mm — l 



r=l s=l 
M-1 



q=0 



r=l q=0 
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that by Eq. (57) becomes 

g=0 s=l 
M-1 

= Y.^r,ih+i){M -l,^)\xr,k), niM <k<{mM-i-l). (58) 

r=l 

The vectors \xr,k) present on the rhs of (58) are hnearly independent for k = 
0, . . . , (M — 1) and the matrix ^V(M — 1,0^ is non singular. Thus, the vectors 
on the Ihs are also linearly independent ioT < h < {M — 1) and rriM k < 
(mjv/_i — 1). The vectors within the sum present on the Ihs of (58) were already 
shown to be linearly independent because they are characterized by q values ranging 
in [0, {rriM — 1)]. One concludes that the vectors \h, k) with < /i < (M — 1) and 
^ k < {rriM-i — 1) are linearly independent. In this way, we have proven 
that the vectors \h,k) with < k < {rriM-i — 1) are linearly independent if < 
h < ifJ'k+i — 1)5 the integers /Xfc being defined by Eq. (39). If rriM-i = rnM, the 
inequality niM < k < rriM-i is never verified and it must be substituted with 
niM^i < k < mM~2 provided mM~2 > '^a/-i(= ttT'm)- In this way, step by step, by 
the procedure just reported one shows that the vectors linearly independent are 

S(a*) = (^\h,k) /i = 0,...,(/ifc-l), A; = 0,...,(mi-l)| 

\h,k) k = 0,...,{mh-l), /i = 0,...,(M- 1)|. (59) 

The corresponding set of {h, k) is nothing else that set X, specified by Eq. (39) 
and shifted by (-1,-1). Hence it is fully specified by the geometry of the 5^ values, 
once these values are mapped into 1. The set i3(a*) consists of J\f points. It will 
be referred to as the principal basic set of vectors along reciprocal crystallographic 
axis a*: basic because i3(a*) determines a complete basis of 1-L{M), linearly related 
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to that formed by \6i),..,\6ji^) and used to define 7i(A/'), and principal (along a*) 
because B{a*) has the largest extension along a* since, at each step of the procedure, 
we tried to include the largest number of reflections {h, k) lying on the rows parallel 
to a*. Any possible confusion being avoided by the context, the set of reflections 
{h, q) with h a.n k obeying to the constraints specified in Eq.(59) will also be denoted 
by i3(a*) and named principal basic set of reflections. Papers I and II showed the 
existence of less elongated basic sets as well as the procedures for singling them out. 
These procedures are based on an "enlargement" method dictated by the Karle- 
Hauptman zeros found during the basic set search. On the contrary, the procedure 
illustrated above only bases on the geometry of the locations of the scattering 
centres. 

B Generalized Vandermonde determinant 

In Appendix A we showed that the vectors \h,k), with h and k obeying Eq. (59), 
are linearly independent because they form the principal basic set of vectors along 
a*. Then, the associated matrix (V) with its elements defined by Eq. (41) must 
be non-singular. The analytical expression of the determinant of this matrix (V) is 
remarkably simple. To get this expression, we note that the full expression of (V) 
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IS 



(V) 



M M2 'm2 





?M„7i 

«^'72\ 


':i''l7ni 

42''21 

':2''22 


^2 ^22 


. £l„™2-l 

';2''21 
42''22 


Af-1 „ 
«1 <^ 
A/ - 1 n 
^2 '521 

M~l n 
«2 '522 


^2 '27712 




42''2m2 


^2 ''27,12 


42''2m2 


A/-1 n 
«2 ''2„7 



A/ ''Ml 

„'"i-i 



'mi 



A/-1 

?i ''i'l 

AJ - 1 „ 
fl '5?2 



fi-i 

^A/ 'a7 2 



?1 ''l 



A/ - 1 TT 



\^A/ ' 



A/"^A/ A/ 'a/' 



'Mm '•Al'Mr 



1 "12-1 A/-1 Q 



(60) 
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The rows of (V) correspond to (r, s) with s = 1, . . . , and r = 1, . . . , M and 
the columns to (p, g) with g = 0, . . . , (mp+i — 1) and p = 0, . . . , (M — 1). The 
determinant of (V) is a homogeneous polynomial in variables {^} and {77}, because 
each term of det(V) has degree Q with respect to variables {77} and degree P in 
the {^'s. Q and P respectively are 



Q = [0 + 1 + ■ ■ ■ + (mi - 1)] + [0 + 1 + ■ ■ ■ + (m2 - 1)] + ■ ■ ■ + 



M 



+ [0 + 1 



m. 



(61) 



k=l 



and 



P = ■ mi + 1 ■ ma + ■ ■ ■ + (M - 1) ■ (62) 

If, whatever i, rji j = rji i with j I, two rows of (V) are equal and the determinant 
will be equal to zero. Thus, one can write 



det(V) 



M 



n n {vi,j2 - Vi,h 

i=l l<Ji<j2<"ii 



(63) 



Since the total degree of the expression inside the square brackets in (63) is Q, one 
concludes that 71 does not depend on variables {rj} so that 71{{C,}, {i]}) = 7^({,^}), 
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and TZii^}) must be a polynomial of degree P. Evaluate now det(V) by considering 
the (mi X nil) minors contained in the first mi rows of (V). From each column of 
the considered minor one first factorizes .^i to an appropriate power. In this way we 
are left with an (mi x mi) Vandermonde-like matrix with elements equal to rjij. It 
is remarked that, whenever mi > m2, a minor is non singular only if it contains the 
(m2 + 1), (m2 + 2), . . . ,mith columns of (V) and, therefore, each complementary 
minor will not present terms rj^j with 2 < r < M and m2 < J < (mi — 1). 
Assume now that ^i = C,j, j being a particular value such that 1 < j < M. We 
have mi > rrij. We can evaluate det(V) by considering all the minors containing 
the rows (l,r) with 1 < r < mi and the rows {j,s) with 1 < s < rrij. Having 
assumed that ,^i = ^j, each column of one of these minors factorizes .^i^ with p in 
[0, . . . , (M — 1)] depending on the considered column. Thus, we are left with an 
(mi + rrij) x (mi + rrij) Vandermonde matrix whose elements are rji p"^ or rjj^g^ with 
I < p < mi, 1, < g < rrij, while r and s belong to [0, . . . , (mi — 1)]. The rank of 
this matrix at most is equal to mi and, therefore, its determinant will have a zero 
of order rrij = min(mi,mj). Hence, det(V) evaluated by this procedure will have a 
zero, related to the fact that ^i = ^j, of order rnj. Assume now that ^2 = for a 
particular j such that 2 < j < M. Before repeating the reasoning made in the case 
^1 = C,j, we imagine of having developed det(V) with respect to the minors contained 
in the first mi rows. As noted above, each complementary minor will no longer 
contain the columns with exponents m2, . . . , (mi — 1). Each of these complementary 
minor can be developed by considering its (m2 + rrij) x (m2 + rrij) minors contained 
in the rows presenting the factors ^2 and C,j- By the same reasoning made above 
for the case ^1 = ^j, one concludes that each (m2 + rrij) x (m2 + rrij) minor has a 
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rank at most equal to m2 and, therefore det(V) must have a zero of order at least 
equal to rrij = min(m2, rrij) when ^2 = 0- concludes that the zero of det(V) is 
at least of the order min(m,mj) when = C,j and one can write that 

mm(m;2,"i(^)" 



det(V) 



n 6. 

.i<;i<«2<Af 



^i({e},w). 



(64) 



The degree of the ^-polynomial inside square brackets is equal to P, so that TZi is 
a polynomial of the only variables {77}. Thus, combining Eq.(63) with (64), one 
finds that 



det(V) = 7^o 



n i^h- 

l<h<l2<M 



min(m;2,mjj' 



M 



n n {vi,j2 - Vi,n) 

.i=l l<ji<j2<mi 



(65) 



where TZq is a simple constant, eventually dependent on the dimensionality of (V). 
Comparing the "diagonal" term Ufiii^/'^)'^' IlfiiUT=iVij^"^ resuhing from the 
calculation of det(V), starting from the explicit expression of (V), with the cor- 
responding term obtained developing the products present in (65), one finds that 
TZq = 1. Thus, the determinant of matrix (V) is 



det(V) 



n 

l<ll<l2<M 



M 



n n {'ni,j2 - Vi,h) 

■«=1 l<ii<i2<mi 



(66) 



By this result it is possible to get the algebraic expression of the determinant of the 
Karle-Hauptman matrix ^P[i3(a*)]^ associated to the principal basic set of vectors 
B{a*). To this aim, we denote the vectors |h) G B{a*) as jh^) with i = 1, . . . ,Af. 
The element of (v[B{a*)]^ is defined as V,^i[B{a*)] = /h,_hr Then, by 

Eq. (18), one finds that 



-i27r(5j-hj 



, z,i=l,...,Af. (67) 
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These can be written, using Eq.s (26) and (27), as 

(p[S(a*)]) = (Vt) . ■ (v) (68) 

It follows that the determinant of the matrix on the Ihs is the product of the 
determinants of the matrices present in the rhs. Hence, by Eq. (66), the determinant 
of the KH matrix associated to B{a*) is 













2 min(mj,mi) 




- M 


2 








n 






n n \Vi/-ni,j 






ij=i J 




-^<^<3<M 








.i=l l<j<£<mi 





This expression applies both to X-ray and to neutrons. Its value is certainly different 
from zero. It is striclty positive in the first case while, in the second case, its sign 
depends on the sign of the first factor related to the product of the charges of the 
J\f scattering centres. 
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